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We review understanding of kinetics of fluid phase separation in various space 
dimensions. Morphological differences, percolating or disconnected, based on over¬ 
all composition in a binary liquid or density in a vapor-liquid system, have been 
pointed out. Depending upon the morphology, various possible mechanisms and 
corresponding theoretical predictions for domain growth are discussed. On com¬ 
putational front, useful models and simulation methodologies have been presented. 
Theoretically predicted growth laws have been tested via molecular dynamics simu¬ 
lations of vapor-liquid transitions. In case of disconnected structure, the mechanism 
has been confirmed directly. This is a brief review on the topic for a special issue on 
coarsening dynamics, expected to appear in Comptes Rendus Physique. 
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I. INTRODUCTION 


Topics related to phase transitions received significant attention over many decades [MU. 
The objective of this review is to discuss developments in the understanding of kinetics of 
phase transitions [ffii with emphasis on phase separating fluid systems Evolution 

of a system from one equilibrium phase to another is a nonequilibrium phenomenon and hap¬ 
pens via nucleation [l]-[6j] and growth of domains. E.g., when a homogeneous binary 

mixture (A + B), liquid or solid, is quenched inside the miscibility gap, it becomes unstable to 
fluctuations and moves towards the coexisting phase-separated state via formation and growth 
of domains rich in A and B particles. For a vapor-liquid transition, the approach to the new 
equilibrium occurs via formation of particle rich and particle poor regions. 

Typically, the growth of domains is a scaling phenomena [9|], i.e., there is self-similarity of 
patterns at two different times, despite a change in length scale £(t), the average size of domains 
at time t. This is captured in the scaling properties [9| 


C{r,t) = C[r/e(t)], 


( 1 ) 


S(k,t) = £(t) d S[k£(t)], 


( 2 ) 


p(£d,t) (3) 

of the two-point equal time correlation function (C), structure factor (5) and domain size 
distribution function (P). In Eqs. (jTH3]) , r is the distance between two points (rj and fj ), k 
is the wave_ vector, is the size of a domain and d is the system dimensionality. There C(x), 
S(y) and P(z) are master functions, independent of time. The above mentioned correlation 
function is calculated as j9[ 

C(r, t) =< V>( U , t) > - < t ) >< t ) >, (4) 

where -0 is the relevant order-parameter for the transition. S(k,t ) is the Fourier transform of 
C(r,t). The scalar notations r and k imply spherical isotropy and are applicable for systems 
without any bias such that the structures are isotropic in statistical sense. The angular brackets 
in Eq.(j4|) are for statistical averaging. 

Typically, £(t) grows with time in power-law fashion as [9 ? ] 

m ~ r. (5) 

The growth exponent a depends upon d, number of order-parameter components, morphology, 
hydrodynamic effects and conservation of order-parameter {§]. In out of equilibrium systems, 
as is clear from Eq. (IT) , the order parameter ^ is a function of space and time. In a mag¬ 
netic system this can be identified with the local magnetization, in a binary mixture with the 
local difference between the concentration of the two species, for a vapor-liquid system with 
the density. Depending upon the type of transition, the total order-parameter (the local value 
integrated over the whole system) may or may not remain same at all times. E.g., in a param¬ 
agnetic to ferromagnetic transition, where, starting from a zero net magnetization, the system 
gets spontaneously magnetized, the order-parameter is a nonconserved quantity. On the other 
hand, for all types of phase separating systems, assuming that there is no chemical conversion 
or flow of material between the system and a reservoir, this is a conserved quantity. In this 
paper we will deal with conserved scalar order parameter. 
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Major fraction of the literature in kinetics of phase transitions is related to the behavior of 
C(r,t) and understanding of the exponent a. However, there exist other interesting aspects, 


e.g., aging 
quantities, e.g 


and persistence 51]. In aging phenomena, typically one studies the two-time 
the order-parameter autocorrelation function |50| 


C W (M„) =< >< >, 


( 6 ) 


where t w and t (> t w ) are respectively the waiting and observation times. Even though the time 
translation invariance is not obeyed by C age (t,t w ) in out of equilibrium systems, this quantity 
is expected to exhibit scaling j50| with respect to t/t w . An objective in studies related to aging 
is to obtain scaling functions for C age or for other two-time quantities of relevance, in different 
types of phase transitions. Persistence, on the other hand, is related to the time dependence of 
the fraction of unaffected local order parameter. 

Each of these aspects are better understood in nonconserved order-parameter case [9]. For 
conserved order-parameter, significant progress has been made only with respect to the growth 


exponent a , for solid binary mixtures [12-|l8|. For the latter, diffusion (via evaporation and 
condensation) is the only transport mechanism and the growth is characterized by a single 
exponent. In fluids, however, no single exponent describes the entire growth process [9j. This is 
due to the faster transport of material, particularly due to the influence of hydrodynamics, at 
late times. Effects of such fast transport may be manifested in different ways depending upon 
system dimensionality and domain pattern. Thus, kinetics of phase separation in fluids is a 
richer and challenging area. In spite of that, reasonable progress has been made (20l425l. f28M49| . 
The aim of the article is to review some of these works in brief. 

The rest of the article is organized as follows. In section II we discuss the theoretical progress. 
Section III is devoted to the discussion of models and computational methods. Representative 
computational results are presented in section IV. Finally, section V concludes the paper with 
a brief summary and outlook. 


II. THEORIES 


In a solid binary mixture, the rate of change of domain size is associated with the chemical 
potential gradient as [l3j dijdt ~ |V/i|. Using the dimensionality of p, as q/T, 7 being the 
interfacial tension and assuming that a gradient exists over the length of the domain size, one 
obtains [13] 


d£ 

dt 


1 

£ 2 ' 


( 7 ) 


Solution of Eq.JTJ) provides a = 1/3, referred to as the Lifshitz-Slyozov (LS) law 12] and is 


understood to be valid for any domain growth occurring via simple diffusive mechanism. 

Early part of fluid phase separation is also dominated by diffusion and thus, is expected 
to have a = 1/3. At late times, hydrodynamic effects become more important, for percolating 
domain morphology. As already mentioned, while for phase separation in solid mixtures system 
dimensionality and domain morphology do not play important role 18], these become crucial 
for selection of mechanism in fluids. I 11 some cases, different mechanisms may give rise to same 
value of the growth exponent. Thus, very direct computational or experimental probes are 
necessary to validate a particular one. Below we will discuss some relevant mechanisms briefly. 

Typically, phase separating systems exhibit interconnected or droplet morphology, depend¬ 
ing upon the overall composition or density of the system. In FigjTJwe show a vapor-liquid co¬ 
existence curve in temperature (T) versus density (p) plane. There the left branch corresponds 
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FIG. 1. Schematic diagram for vapor-liquid coexistence in T vs p plane. The critical point is marked 
by a cross. The left and right branches represent respectively the vapor and the liquid densities. From 
Molecular simulation. 


to the vapor-phase density and the right one to the liquid-phase density. For a binary mixture 
phase separation, the appropriate abscissa variable is the concentration of one of the species. 
For density or composition very close to the coexistence curve (vapor branch in vapor-liquid 
transition), droplet like pattern is visible, from the very beginning. In this under-saturation 
case, the system requires fluctuations involving large length scales to nucleate stable droplets 
that can grow. Such events, depending upon the proximity to the coexistence curve, may be 
rare and the phase separation can be delayed. On the other hand, far from the coexistence 
curve, the system falls unstable to infinitesimal fluctuations and the corresponding phase sepa¬ 
ration process is referred to as the spinodal decomposition. Typically, in this case, the domain 
morphology is percolating. Here we caution the reader that at late time droplet structures 
can emerge from percolating morphology as well [36j, if the system is away from the critical 
composition or density. In fluids this is expected to happen faster. 

Let us first consider the case of percolating domain morphology. Because of experimental 
realizability, we will first take up the d = 3 case. The hydrodynamic regime is divided into 
two subregime s (p , viz., viscous and inertial. These we discuss below, following Siggia [ 22 f| and 
Furukawa [23|, (24|], using the Navier-Stokes (NS) equation 

p—v — uV 2 v = — Vp. ( 8 ) 

In Eq.(JHJ), p is the mass density, v is the fluid velocity, v is the kinematic viscosity, p is the 
pressure and D/Dt = d/dt + (F.V). In Dv/Dt , both the terms stand for acceleration, as can 
be checked dimensionally. The second one is related to convection where the acceleration is 
generated by space dependent velocity. Eq. (jHJ) is essentially related to balancing of force coming 
from pressure with that of frictions. On the left hand side, the first term can be treated as 
inertial friction and the second one is related to dissipative friction originating from viscous 
drag. Via dimensional substitutions Dv/Dt = v/t, V 2 v = —v/£ 2 and considering that the 
pressure gradient is obtainable from interfacial tension 7 , one writes, after multiplying both 
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sides by vt, 


2 P u 7 

P v + “ T v = 4 - 


( 9 ) 


The second term on the left hand side of Eq.()9l) is related to the inverse of time dependent 
Reynolds number and is important at earlier part of the hydrodynamic regime, providing viscous 
growth. At later time, the first term dominates, giving rise to the inertial hydrodynamic growth. 
Thus, in the first part of hydrodynamic growth it is equivalent to equating the interfacial free 
energy density with the viscous stress and in the second part, the same with the kinetic energy 
density [22|. With the understanding that there exist unique time and length scales in the 
problem, one identifies v with di/dt . Then, in the viscous regime one obtains 


and in the inertial regime 


di 

dt 


M 

dt 


pv 



( 10 ) 


(ii) 


Solutions of Eqs. (ITO]) and (ITT]) provide a = 1 and 2/3, respectively. 

Recall that our discussion is related to interconnected domain morphology, having tube-like 
structure in d = 3. The fast advective transport of material through these tubes happens due 
to pressure gradient calculable from interfacial tension and domain size. In d — 2, analogue 
of a tube is strip. San Miguel et al. 25[ argued that fluctuations along the linear boundaries 
of these strips significantly increase the interfacial free energy due to increase in curvature. 
This mechanism thus lacks instability. According to these authors, instead of a crossover from 
a = 1/3 to 1, as in d — 3, the 2 — d fluid will encounter a crossover from 1/3 to 1/2, the 
latter coming from an interface diffusion mechanism [25[. There a crossover from a = 1/3 to 
1 will require large fluctuations, may become possible only at very high temperatures. In the 
interface diffusion mechanism, break-up of strips into droplets is expected 125]. In that case, 
droplet diffusion and collision mechanism of Binder and Stauffer (BS) [20H22J, in general valid 
for off-critical fluid quenches, is expected to be operative. This we discuss below. Note that, 
in case of droplet morphology, the growth description cannot be provided via NS equation [§]. 

As opposed to solid binary mixtures, droplets in fluids are significantly mobile. This 
gives rise to, in addition to the standard evaporation-condensation process, another growth 
mechanism involving diffusion and collision of droplets. These are inelastic collisions, following 
which the colliding partners merge, forming a bigger droplet and reducing the droplet density 
For the decay of n, one can write j22j 


[n 


dn 

dt 


= KD e Zn 2 , 


( 12 ) 


where K is a constant and Dp is the diffusion constant for a droplet of size t. Following Stokes- 
Einstein-Sutherland relation 52], Dpt can be treated as a constant. Using n ~ l/t d , one obtains 


22 


dt 


dt t d 


oc 


-l 


(13) 


Solution of Eq. (1 13 1) provides a = 1/d. Thus, in d — 2, BS value of a is same as that for the 
interface diffusion mechanism of San Miguel et al. On the other hand, in d — 3, BS mechanism 
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provides a value of a coinciding with the LS value. It is expected that the amplitudes will be 
different for different mechanisms that provide same exponents. E.g., in d — 3 LS mechanism 
should have a smaller amplitude than the BS one. Direct verification of the latter is possible via 
computer simulations 45]. Tanaka 2^-2^] argued that in high droplet density situation, motion 


of the droplets will not be random due to inter-droplet interaction. However, incorporation of 
such fact merely changes the growth amplitude. Presence of such interaction can be verified by 
calculating the mean-squared-displacements (MSD) [52] of the droplets. 

Pictures presented above are primarily meant for kinetics of incompressible liquid-liquid 
phase separation. Nevertheless, they may apply for vapor-liquid phase separation to a good 
extent. This, however, requires verification. 


III. MODELS AND METHODS 


Kinetics of phase separation in solid binary mixtures at atomistic level have been tradition¬ 
ally studied via Monte Carlo (MC) simulations [53] of the nearest neighbor Ising model 


H = 


-J £ S.Sy, 

<ij> 


J > 0; Si = ±1, 


(14) 


on regular lattice. Here, an up spin (S l = +1) corresponds to an A particle and a down spin 

a = — 1) to a B particle. In this model, the kinetics is introduced via the Kawasaki exchange 
trial moves in which positions of two nearest neighbor particles are interchanged. These 
moves are accepted according to the standard Metropolis algorithm [53j. At the coarse-grained 
level, one solves the Cahn-Hilliard (CH) equation [9] 


dip{r,t) 

dt 


= -V 2 


ip(r, t ) + V 2 ^(r, t) — ip 3 (r, t) 


(15) 


where the order-parameter ip can be thought of being obtained from the coarse-graining of the 
Ising spins over a length scale of the order of the equilibrium correlation length £. Thus, the 
CH equation provides an advantage of accessing large effective length scale within practically 
accessible computation time. Eq. (ll5j) can one way be obtained from the Ising model via a 


master equation approach with Kawasaki kinetics as an input [54], [55]. Despite the latter being 


a direct and rigorous approach, in this article we will present a phenomenological derivation. 
Eq. (fT5]l is solved on regular lattice, typically via Eular discretization technique [55j. This model, 
according to the nomenclature of Hohenberg and Halperin []56j, is referred to as the model B. 
In fluids, at coarse-grained level the kinetics is typically studied via model H 


56 


which is 

a combination of the CH and the NS equations, the latter taking care of the fluid flow. To 
describe this model, derivation of the CH equation (9] will be of help. Recalling that we are 
dealing with conserved order-parameter, one writes the continuity equation 


dip 

~di 


= -v ■ J c , 


(16) 


where the current J c comes from the gradient of the chemical potential — V/i. The chemical 
potential can be obtained from the functional derivative as 


h = 


8F[rp(r,t)\ 


(17) 
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the free energy functional being the Ginzburg-Landau one 

F = k B T / df[-aip 2 + bip‘ + e(W>) 2 ], 


(18) 


where the positive coefficients a, b, and c are temperature dependent. Often these coefficients 
are scaled in such a way that the final form is parameter free, as Eq. (fT5]b 

The model H equations read j9;], in the incompressible limit (V • v — 0), as 




(19) 


( 20 ) 


Dv 2 _ 

p-£j[ ~ v = — wVp. 

In addition to the pure CH and NS equations, there are additional terms in Eqs. (fl9|) and 
coming from the coupling between the velocity and order-parameter fields. This is justifiable, 
e.g., the current J c in the CH equation is expected to be influenced by the velocity held, 
becoming 


J c = -V/j - tpv, 


( 21 ) 


and the chemical potential gradient should provide a driving force to affect the velocity held 
in the NS equation. Again, these equations can be solved on regular lattice. This model is 
very much phenomenological and a microscopic derivation like the CH equation [HU] is still 


lacking. Such an objective can possibly be achieved 57] via the construction of free energy 


functionals by taking density and velocity distributions as inputs from molecular dynamics 
(MD) simulations at various coarse-grained levels, by appropriately adjusting the effective inter¬ 
particle potential at successive steps. A multi-scale method like this may seem obvious but has 
never been demonstrated. Success of the method may prove useful in justifying similar (elegant) 
phenomenological models in more complex problems like active matter |58| ]. 

The incorporation of the coupled terms in Eqs. (TT9|) and (1201) . in a sense, guarantee the 
simultaneous growth in the velocity and density helds. While such facts have been reported 


from Lattice Boltzmann simulations [10], ]34J, are lacking in MD simulations even though the 
latter is more direct and accurate. Thus, a microscopic justification of the model, particularly 
understanding of the importance of additional higher order terms, is needed. 

As can be guessed from the above discussion, at atomistic level, kinetics of fluid phase 
separation is studied via MD simulations [56j, [59|, 60] in which hydrodynamics can be easily 
implemented. For this purpose, the interatomic interactions are popularly incorporated via the 
well known Lennard-Jones (LJ) potential 


u(r) = 4e 


O-N 12 

r, 


r 


( 22 ) 


where r (= \r\ — fj |) is the distance between two particles at r % and rj, e is the interaction 
strength and a is the interparticle diameter. In a binary (or multicomponent) mixture, £ and 
a can be chosen to be different for different combinations of species. 

In standard MD simulations, equations of motion are solved in discrete time, typically via 
Verlet velocity algorithm [56=, 59]. Conservation laws of hydrodynamics are perfectly satisfied in 


microcanonical ensemble using which various fluid transport properties are calculated in equi¬ 
librium. However, for kinetics of phase separation, following a temperature quench, in course of 
the system’s evolution, the potential energy decreases, for energy driven phase transitions. In 












a constant energy ensemble, thus, the kinetic energy, i.e., the temperature increases, destroy¬ 
ing the objective of the study. Therefore, simulations in canonical ensemble, with appropriate 
hydrodynamics preserving thermostat, become essential. There are a number of good methods 
for this purpose, e.g., dissipative particle dynamics [61 63 
Lowe-Andersen thermostat 6 f 


multiparticle collision dynamics 
Nose-Hoover thermostat §39], etc. In this article, pre¬ 
sented results were obtained via application of the Nose-Hoover thermostat (NHT) which will 
be discussed soon. 


IV. SIMULATION RESULTS 



U(r < r c ) 
U(r > r c ) 


du(r) 

u(r) - u(r c ) - (r - r c )—— 

CLT r=r, 

o. 


(23) 


The cut-off radius r c (= 2.5 a) is introduced for a faster computation. Note that, in critical 
phenomena 0,0, LJ potential being a short-range one, this does not alter the universality 
class. Due to the cut-off and shifting of the potential to 0 at r = r c , a discontinuity in the 
force is created which may cause non-smooth behavior in energy and thus may be problematic 
for energy and momentum conservations, necessary for preservation of hydrodynamics. This 
problem is corrected by incorporating the last term in the first part of Eq. (1231) . 

The phase diagram, the primary requirement before studying the kinetics of phase sepa¬ 
ration, for atomistic models can be obtained using both MC and MD simulations. As is well 
known, equilibrium phase behavior and thermodynamic properties are insensitive to ensemble 
and technicality related to transport mechanism. Thus, MC simulations with smart ensem¬ 
bles are advantageous, compared to MD. For vapor-liquid phase transition, typically one uses 
grandcanonical and “Gibbs” ensemble methods [5.3]. With grandcanonical ensemble, other 
thermodynamic properties, e.g., compressibility, can be easily obtained to study its critical 
singularity. 

As mentioned, an NHT was applied to study dynamics via MD. There one solves the 
deterministic equations of motion 59| 


mri = pi , 


Pi = - 7 ^ - -Pi, 

OT i 


(24) 

(25) 


£(p?/m) - 3Nk B T 


Q 


( 26 ) 
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In Eqs.dMHm is the particle mass (set equal for all), 5 is a time dependent drag that 
adjusts its value depending upon the drift of temperature from the assigned value, N is the 
number of particles, ks is the Boltzmann constant and Q is the coupling strength between 
the system and the thermostat. Essentially, in this scheme one works with a microcanonical 
ensemble with a modified Hamiltonian. Results thus obtained are equivalent to those from the 
canonical ensemble with the original Hamiltonian 59, 162 


There exist better hydrodynamics preserving thermostats. But for the present purpose, the 
NHT appears adequate. The results obtained using this thermostat have been tested against 
other methods. For example, in d = 2, we have compared j4l| them with the Lowe-Andersen 
thermostat (LAT). The LAT is an improvement over the basic Andersen thermostat (AT). 
In the AT randomly chosen particles are assigned new velocities (with Maxwell distribution) 
to keep the temperature at a desired value and thus stochastic in nature, like MC methods. 
In the LAT, while assigning new velocities, conservation of local momentum is appropriately 
taken care of. Further, inside the coexistence region, transport properties of droplets (from 
equilibrium configuration), calculated via NHT, are found to be in good agreement with the 
calculations using microcanonical ensemble 46, 47] . Also, in a binary fluid, critical behavior 


of shear viscosity was nicely reproduced via NHT 6 f|. In addition, from experience we feel 
that the NHT provides a superior control over temperature, compared to few other better 
hydrodynamics preserving thermostats, particularly in out-of-equilibrium situation. 

We will present results 42-49| in both d — 2 and 3, from MD simulations in periodic 
square or cubic boxes of linear dimension L (in units of cr). All results were obtained after 
averaging over adequately large numbers of independent initial configurations. In the solutions 
of the equations of motion we have used time discretization At = 0.005f 0 , t 0 [= (mo 2 / e) 1 ^ 2 ) 
being an LJ time unit. Note that the critical temperature and the critical number density for 
this model in d — 3 respectively are ~ 0.9 e/ks and ~ 0.3. All results in d = 3 correspond to 
quenches from T = oo to 0.6 e/ks- As expected, the value of T c in d = 2 is lower (~ 0.5 e/ks)- 
The quench temperature there is 0.25 e/ks- 

In Figj2]we present snapshots from d = 3 (unless otherwise mentioned, results are from 
this dimension) for different overall densities . It is clearly seen that as one approaches the 
vapor branch of the coexistence curve, the morphology becomes more droplet like. Given that 
we have access to both percolating and droplet structures, a wide variety of mechanisms can 
be checked. 

In Figj3]we show a scaling plot of the correlation functions for droplet morphology, using 
data from different times. For this calculation, the order parameter -0 was given a value +1 if 
the density at a point is higher than the overall density p and —1 otherwise. Nice collapse of 
data, versus r/£(t), confirms self-similarity of structure. Note here that we have obtained £(t) 
from different methods, from the lengths at which C(r,t ) decays to a particular value, from 
the first moment of P(£d,t), as well as by identifying domain boundaries and directly counting 
the number of particles inside them. For the last method, cube (square in d = 2) root of the 
average number of particles inside the liquid regions will provide £(t). For the sake of brevity 
we present results only from the correlation functions. 

In Figj4]we show plots of £(t) vs t, for various values of p. It appears that closer to the critical 
density there is a linear growth, after a brief slow regime, before hitting the finite-size effects. 
The slower part at the very beginning can be attributed to the LS behavior [12] and the linear 
one to the viscous hydrodynamic one 22-24], Due to the demanding nature of MD simulations, 


to the best of our knowledge, a crossover from viscous to inertial hydrodynamic regime is not yet 
appropriately observed, using this method. However, for intermediate densities, we observed a 
f 2 / 3 behavior even before a linear growth appears. This has been checked [47[ via appropriate 
f ini te-size scaling analysis j67j which we do not present here. A possible reason for observing 
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FIG. 2. Evolution snapshots for various overall densities p , obtained from molecular dynamics simu¬ 
lations of the single component Lennard-Jones model in d = 3. In all the cases the time, temperature 
and the system size were fixed to t = 150, T = 0.6 and L = 80. From J. Chem. Phys. 139 , 044911 
(2013). 



FIG. 3. Scaling plot of C(r,t ) as a function of r/t(t ), for p = 0.05. Data from few different times are 
included. The results correspond to droplet morphology in d = 3. From Soft Matter 9 , 4178 (2013). 


the exponent 2/3 during early period of the hydrodynamic growth can be the following. In the 
inertial regime, due to the large mass of domains, there may be competition between growth and 
break-up. In intermediate density regime, where interconnectedness of the domain morphology 
is not very robust, the break-up becomes easier. For densities very close to the coexistence 
vapor density, certainly the growth is much slower. This we discuss in details below (45-471 . 

In Fig|5]we take a look at the results for p = 0.05 on a double-log scale. The flat behavior 
of the data at the beginning is due to delayed nucleation. A very fast rise after this signals the 
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FIG. 4. Plots of £{t) versus t, for different values of p, in d = 3. The continuous straight line 
corresponds to viscous hydrodynamics growth. From J. Chem. Phys. 139 , 044911 (2013). 



FIG. 5. Log-log plot of £ versus t, for p = 0.05, that corresponds to droplet morphology in d = 3. The 
solid lines correspond to power-laws, exponents being mentioned on the figure. From Phys. Rev. E 
85 , 050602(R) (2012). 


onset of nucleation, following which the data are consistent with a t 1 / 3 
already, in d = 3 this behavior can be due to LS 


13 as well as BS [20-22 


;rowth. As mentioned 
mechanism. To 


select from the multiple possibilities we do the following exercise. 

In Figj6]we show plots of dn/dt versus n, for two different low values of p, both providing 
droplet structure, on a double-log scale. Technical details on the counting of droplets by 
identifying them can be found elsewhere (45-4^1 • Linear look in both the cases confirm power- 
law behavior. The solid lines there are proportional to n 2 with which the simulation data are 
consistent. This verifies the starting equation [see Eq. (jl2p ] for deriving the BS 20-22] growth 
law. Thus, the possibility of a LS 12[ mechanism is ruled out. We have also checked that 
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FIG. 6. Plots of dn/dt versus to, n being the droplet density. Results from two values of p, in d = 3, 
are included. The solid lines are proportional to to 2 . For p = 0.08 a multiplicative factor 5 was used 
to separate the data sets from each other. From J. Chern. Phys. 139 , 044911 (2013). 



FIG. 7. Mean squared displacements for droplets in equilibrium and nonequilibrium situations are 
plotted versus time. The droplets are of similar size. The quadratic and linear behavior, corresponding 
to ballistic and diffusive regimes, are appropriately marked by different lines. The results corresponds 
to d = 3. From Soft Matter 9 , 4178 (2013). 


between two collisions, sizes of droplets remain same, within minor fluctuations. This further 
discards the importance of LS mechanism in the present context. 

Next, we come to the point of inter-droplet interaction. To understand it we have calculated 


the MSD [52| of individual droplets [H, 0 


as 


MSD =< (R CM {t) - Rcm(0)) 2 >, 


(27) 


where Rqm (t) is the location of the centre of mass of the droplet under consideration at time t. 
In Fig. [7] we show MSD for a typical droplet during nonequilibrium evolution, as a function of 







13 



FIG. 8. Plot of i vs t, obtained from the simulations of LJ fluid in d = 2. The solid line has i 1 / 2 
behavior. The results correspond to droplet morphology for overall density p = 0.02. 


time. Note that the duration over which such data can be presented is dictated by the average 
collision interval. This is larger for low droplet density and increases with the evolution of 
the system. In the same graph we showed a corresponding plot for a droplet in equilibrium 
situation, of approximately the same size as the nonequilibrium case. The difference between 
the two cases is significant. On a double-log scale, at late time, the equilibrium droplet exhibits 
linear behavior, confirming Brownian motion. The super-linear behavior in the nonequilibrium 


26-28 


case is due to inter-droplet interaction, as stated by Tanaka 

Finally, in FigJH] we show i versus t plot in d — 2 for off-critical composition 


Here 


clearly an exponent a = 1/2 is visible. From the calculation of dn/dt , in this dimension also we 
have confirmed that the mechanism is BS. For the sake of brevity we do not present these here. 
For higher overall density, the interface diffusion mechanism in this dimension was previously 
observed 


V. CONCLUSION 


We have provided a brief review on kinetics of fluid phase separation. A general discussion 
on various theoretical pictures, based on the system dimensionality and domain morphology, is 
provided. Methodologies related to coarse-grained and atomistic models are discussed. Partic¬ 
ular emphasis was on the molecular dynamics simulation methods. 

Molecular dynamics results for vapor-liquid transitions are presented (421- 1471 Sj in d — 2 
and 3. These include percolating as well as droplet structures. The domain growths from these 
simulations are observed to be consistent with theoretical predictions, despite the fact that 
most of these predictions are related to phase separation in binary fluids. For brevity, we have 
avoided results on binary fluid mixtures. 

Recently, important results have been obtained with respect to aging in fluid phase sep¬ 
aration, for binary fluids as well as vapor-liquid systems, via molecular dynamics simulations 
68 ;, .69}. These results are understood via simple scaling arguments. In addition, effects of disor¬ 
der in kinetics of fluid phase separation have been looked at [4l[ and compared with disordered 
Ising systems. Lack of space, however, restricts us from including these results. 
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Further, there have been significant activities in the area of kinetics in confined geometry 


20-480]- In computational front, good number of reports exist in fluid phase separation 74 


0 7§, |80| as well. However, the understanding of these results are not as complete as pure 
two- and three- dimensional systems. 

Even though there exists significant agreement between theories and simulations of simple 
models, the situation is not as satisfactory as far as experiments are concerned. Discrepancy 
between theoretical predictions and experimental observations can be attributed to the presence 
of impurities in real systems as well as to over-simplified pictures in theoretical calculations. 
On the other hand, it is possible to do simulation study of more realistic models |4lL 81 1 to 


obtain better agreement with experiments, for both diffusive and hydrodynamic coarsening, 
particularly when improved computational resources are available. 
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